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Abstract. We derive oscillatory signals in correlation functions in two-field open inflation by 
means of the in-in formalism; such signatures are caused by resonances between oscillations 
in the tunnelling field and fluctuations in the inflaton during the curvature dominated, in- 
termediate and subsequent inflationary regime. While amplitudes are model-dependent, we 
find distinct oscillations in the power and bi-spectrum that can act as a direct probe of the 
curvature dominated phase and thus, indirectly, strengthen the claim of the string landscape 
if they were observed. We comment on the prospects of detecting these tell-tale signs in 
current experiments, which is challenging, but not impossible. 

At the technical level, we pay special attention to the applicability conditions for trun- 
cation fluctuations to the light (inflaton) field and derive upper limits on the oscillation 
amplitude of the heavy field. A violation of these bounds requires a multi-field analysis at 
the perturbed level. 
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1 Introduction 

The existence of a landscape in string theory with a plethora of metastable vacua suggests 
the presence of eternal inflation, yet our universe must have undergone around sixty e-folds of 
subsequent slow-roll inflation to solve the problems of the big bang and lead to the observed 
fluctuations in the cosmic-microwave background radiation. These two regimes of inflation 
are connected via quantum tunnelling. If the transition can be described by a Coleman- 
De Luccia instanton [1], the universe in the nucleated bubble appears as an infinite, open 
universe. Spatial curvature is subsequently reduced below current observational bounds by 
slow-roll inflation. Hence, this model of open inflation [2-6] seems to be a generic occurrence 
in string theory, providing (in principle) a theory of the initial conditions for the slow- 
roll phase at the background and the perturbed level. Two observational tell-tale signs 
of open inflation are well known: spatial curvature, see e.g. [7] for a recent analysis, and 
super-curvature fluctuations of varying origin [8-11] (see also [12-16]) that affect correlation 
functions on the largest scales [3]. 

In this paper, we investigate a third indicator that is generically present if the tunnelling 
field a is not identical to the inflaton as in [4, 5], leading to quasi-open inflation [6]: since a is 
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usually displaced from the inflationary valley after tunnelling, as in the concrete case study 
of [17], or because it got excited subsequently, it exhibits oscillations during the curvature 
dominated phase that result in additional oscillatory signatures in correlation functions, 
whose frequency in k-space, u}{k), acts as a direct probe of the background scale factor and 
thus the curvature dominated regime. 

The use of signals induced by oscillating heavy fields to identify the scale factor's evo- 
lution was first advocated by Chen [18-21], particularly to differentiate slow-roll inflation 
from proposed alternatives, such as bouncing cosmologies. We extend Chen's results to open 
inflation, for which ome of the approximations made in [19, 20] are not applicable. 

After introducing a simple two-field open inflation model [17] as the background of 
interest in Sec. 2, we turn to perturbations: following the philosophy of [19, 20], we truncate 
perturbations to fluctuations in the inflaton, and treat oscillations in cr as a small correction 
to the background. We pay special attention to the validity range of such a truncation, 
resulting in strict upper bounds on the allowed oscillation amplitude of a, see Sec. 3.1. These 
bounds are often not stated explicitly in the literature (see however [22, 23]). While some 
case studies in [19, 21] contain amplitudes in excess of our bounds, the qualitative point of 
these studies remains valid. If larger amplitudes are to be considered, the full multi-fleld 
setup needs to be solved at the perturbed level. The effective-fleld theory proposed in [24- 
29] has a different validity range, which does not cover the regime we are interested in (see 
Sec. 3.1 for a brief discussion). 

Given such a truncation, we use the in-in formalism [30], briefly summarized in Sec. 3.3 
(see [31] and [32] for pedagogical reviews), to compute correlation functions in Sec. 3.4: 
resonances between the mode functions and oscillations in background quantities caused 
by a{t) lead to oscillatory signals. Since these resonances concern only perturbations deep 
inside the horizon, we are justified to use the Bunch-Davies state as an approximation to the 
(excited) state after tunnelling. Sec. 3.2. 

Based on the stationary phase approximation, we derive analytic expressions for the 
power-spectrum and the bi-spectrum in three smoothly connected regimes, curvature dom- 
ination (Sec. 4.1), inflation (Sec. 4.2) and an intermediate regime (Sec. 4.3). We recover 
Chen's results during inflation and provide concrete case studies for signals generated during 
the prior two regimes. Sec. 4.1.1 and Sec. 4.3.1. 

We find that while an initial displacement of a from the infiationary valley directly after 
tunnelling leads to distinct oscillations in correlations functions, their amplitude is below 
observational sensitivity in both the power- and bi-spectrum in current experiments (within 
our computational strategy). However, a subsequent displacement of a in the curvature 
dominated regime, which can take place without spoiling the validity of our approximations 
(Sec. 3.1.1), can lead to signals at the percent level in the power-spectrum, if conditions 
are favourable (u is heavy and it is displaced shortly before the intermediate regime). Such 
signals are within reach of observations and should be sought after. The amplitude of the 
bi-spectrum remains small. We discuss our results in Sec. 5. 

2 Two-field Open Infiation 

In open infiation [2, 3, 3-5], the slow-roll regime is preceded by a false vacuum driven phase 
(eternal infiation), which is terminated by a Coleman-De Luccia tunnelling event [1] (see [33] 
for a textbook treatment). From within the nucleated bubble, the Universe can be described 
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as an infinite, open FRW universe ^. Around 60 e-folds of slow-roll inflation are needed to 
reduce curvature below current bounds, = O-OO^T^QQQgg [34] (degeneracies with regard to 
the type of dark energy are lifted by current data [35]). 

False vacua are abundant on the landscape of string theory [36] , rendering open inflation 
an appealing candidate on this higher dimensional field space with around 500 moduli fields. 
After the tunnelling event, the fields do not find themselves generically in another local 
minimum; if the target region happens to be flat enough, a subsequent phase of slow-roll 
inflation can result. 

From a phenomenological point of view, both large- and small-field models are possible; 
the latter ones may be more appealing, since they avoid the rj problem and appear to be 
natural on "random" landscapes, for which slow- roll infiation usually takes place in the 
vicinity of an infiection point, see e.g. [37-42]; a short duration of inflation is more likely 
than longer ones [37, 38] (the probability often scales as A^~^ [42, 43], where is the number 
of e-folds). If this phase lasts less than around sixty e-folds, curvature disrupts subsequent 
structure formation, preventing the genesis of galaxies and thus observers [43]. Hence, based 
on mild anthropic reasoning, we expect curvature to be close to observational bounds [34, 35] , 
leading to an indirect observational test of the landscape idea [7]^. 

In addition to searching for spatial curvature directly, one may observe other tell-tale 
signs of open inflation: slight deviations from a nearly scale invariant power-spectrum on the 
largest scales, which arise due to the presence of super-curvature perturbations [8-16]. These 
deviations trace back to the difference of the state from the Bunch Davies vacuum after the 
tunnelling event [9]. They might help to explain the observed lack of power on the largest 
scales, but a conclusive observation is lacking due to cosmic variance. 

It is also possible that our bubble collided with other ones during the slow-roll regime 
[49], which may lead to giant circular structure in the CMB, see e.g. [50, 51] among others. 
Such peculiar patterns have not been found beyond the 2a level [52]. 

It is possible to realise the one-bubble open inflation scenario with a single scalar field 
[2] responsible for both, the tunnelling event and subsequently inflation. However, functional 
flne-tuning is needed, since slow-roll inflation requires a small fleld mass, whereas tunnelling 
via a Coleman-De Luccia (CdL) instanton hinges on a steep, thin barrier. Further, in light 
of the many fields present on the landscape, it seems more natural to employ at least two 
fields, one responsible for tunnelling and a second (perpendicular) one to instigate slow-roll 
inflation. Naturally, even more flelds could be invoked, for instance to end inflation, generate 
fluctuations via the curvaton mechanism or generate oscillatory signals during the slow-roll 
regime by means of turns in the trajectory [23, 53-58, 60]. (see [24-29] for a series of papers 
developing an effective field theory framework). 

Since we are primarily interested in additional observational handles on open infiation, 

^In the original treatment by Coleman and De Luccia (CdL), the bubble expands almost at the speed of light 
after it has nucleated; the inner edge of its wall traces out hyperbolic surfaces in the surrounding Minkowski 
space (in the bubble interior) which are asymptotic to the light cone at the origin of the bubble. These 
surfaces of constant (negative) spatial curvature are precisely where the field is constant. They subsequently 
correspond to surfaces of constant CMB temperature, which is our definition of constant cosmic time. Thus, 
the whole of the observable (infinite) universe today fits in the space-time interior of the finite bubble. In 
the two-field model that we consider the universe inside the bubble is actually finite [6], which is sometimes 
referred to as quasi-open infiation. 

^This line of reasoning is hampered by the measure problem [44, 45] (see [46] for a review or proposed 
measures). The earliest attempt at making quantitative predictions for ilk in open infiation can be found in 
[47], and for quasi-open infiation [6] in [48]. 
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we focus on a simple two-field model (without turn): a denotes the direction in which tun- 
nelling takes place, and a perpendicular field, (j)^ subsequently drives inflation. A concrete 
computation of the CdL instandon in such a two field model can be found in [17]. 

2.1 The Setup 

Consider two canonically normalized scalar fields with the action 

I' ^ ( Ml, I 1 \ 

S = j d^x^g I -^R - -^g^^d^od^a - -g^''d^ct>dA - W[cj, 0) j . (2.1) 

Given the desired distinction between a and we take a separable potential 

W = V{<j,) + ^m^a^ , (2.2) 

where V{(f)) is assumed to be sufficiently ffat to yield slow-roll inflation (in this study, we 
won't need to specify its exact form for the most part). The metric of an open, FRW universe 
is 

^''^ +r2(d02^sin2 0d(/.2) 



ds'^ = -dr + 



l-kr^ 



(2.3) 



with k = —1 and we use the notation (...)' = d{. . .)/dt. Occasionally, we use conformal 
time, adr = dt, to simplify analytic expressions. 

We treat o" as a small, homogeneous perturbation to the background, so that we can 
define a background Hubble parameter via 

which dominates the Hubble friction in the equation of motion of a, 

a + 3Fbgr<7 + mla w . (2.5) 
To ease notation, we set the reduced Planck mass to unity throughout this article, 

M„, = -^.l. (2.6) 

We note that a ^ after tunnelling [17], but even if it were, it could get displaced due 
to a feature in the inflationary attractor valley, see Sec. 3.1.1 for a brief discussion; we do 
not explicitly incorporate features in our potential, but use this possibility as a motivation 
to treat the initial amplitude (Tq as a free parameter, set either directly after tunnelling, or 
at a later time. The maximal oscillation amplitude is bounded from above, if a truncated 
single-fleld description is to be used at the perturbed level; we elaborate on this subtle issue 
in Sec. 3.1. Any displacement of a leads to damped oscillations and thus a contribution to 
the Hubble parameter 

^' = y + ^bgr' (2-7) 

of 

Pa = l{a' + mla^)^3H^,. (2.8) 
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To treat perturbatively, we need 



« 1, (2.9) 

P0 + Pk 

where we defined 

Pk = — = ^, 2.10 

which leads to an upper bound on do that is, however, weaker than the one we impose later 
on to stay within the regime of a truncated single-field description at the perturbed level, see 
Sec. 3.1.1. By choosing the normalization k = —1, the scale factor has units of 1/mass. 

Since we are primarily interested in the transition from the curvature dominated regime 
to inflation, we sometimes approximate 

^ const (2-11) 

at the background level ^. 

We will need the leading order oscillations of the Hubble parameter, Hosd, to compute 
corrections to correlation functions of the curvature fluctuation C; these are caused by reso- 
nances with mode functions in the in-in formalism, see Sec. 3.4. To this end we separate the 
Hubble parameter into a smoothly varying part, and an oscillatory one 

H = Ho + Hosci, (2.12) 

where Hosd ^ Hq. Since the only oscillatory contribution stems from (7{t), we can identify 
-f^osci by comparing (2.7) with 

^ Hq + IHqHoscI ■ (2.13) 

Since a closed form expression for H„ is not attainable in the transition region to inflation, 
we apply a moving average to (2.8) in order to extract the average and the corresponding 
oscillatory component numerically in this case. Given i/osci) the leading order oscillatory 
contributions to combinations of the slow-roll parameters 



IL -A. 



that we need later on become 



fr2 



(2.14) 



(2.15) 



I f ... 3HH 4i7M 



^We do not suggest that inflation is driven by a cosmological constant, but merely make this approximation 
for certain background quantities. 
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2.2 From Curvature Domination to Inflation 

We can identify three different regimes in tlie evolution, 

1. curvature domination, S> {Pcf>,Pa}, 

2. an intermediate regime (transition to inflation), Pk P<j) ^ Pa , 

3. and inflation, p^f, » {pk,Pa}, 

each of which we can treat separately. Since simple analytic expressions are often attainable 
in conformal time, we shall freely switch between t and r in the following. We put the 
tunnelling event in the far past, rtun 0. 

Since spatial curvature dominates the energy balance initially, the scale factor grows 
exponentially with r. We start our investigation at 



(2.17) 



so that the curvature contribution to the Priedmann equation is unity in natural units. This 
phase ends once the curvature contribution approaches that of the inflaton at 

of) = ^ , (2.18) 

with ^ > 1 (we choose A = 10). As the expansion progresses, energy stored in spatial 

curvature continues to red-shift oc a~^, while p^ ~ const, and p^ oc (a oscillates in a 
quadratic potential). After a brief intermediate regime, inflation commences, during which 
we approximate aoct^ with j9 ^ 1. We define 

(2.19) 



with i? < 1 so that (j) dominates for a > af (we choose B = 1/A = 0.1). We denote the 
instance of equality by Tgq, when = H^, or equivalently 

« = «eq = . (2.20) 

We count the number of e-folds positive from here on, 

/•*end 

N= Hdt, (2.21) 

and assume a V{4>) such that inflation ends once A'^ > 60 and the remaining curvature is 
diluted sufficiently to solve the flatness problem. The length of slow-roll inflation can be 
tuned freely by adjusting 1^(0) appropriately. At the end of inflation, any cr-oscillations have 
died away, since a is assumed to be considerably heavier than the inflation. During inflation, 
the leading order term to the slow-roll parameter is constant and small, e = p~^. 

Given these three regimes, it is simple to derive approximate analytic expressions at 
the background level. A particular combination that enters our subsequent computations is 
defined by 

/(r) = i^ + J-^^, (2.22) 
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and a straightforward computation, entailing approximate solutions to the Friedmann equa- 
tion, yields 



/(r)ss< -tanh(r) for a' 

^ N for a 

\ (i-p)t 




(2.23) 



The to corresponds to the choice of initial time, see Sec. 4.1. 
3 Resonant Signals in the In-In Formalism 

We want to compute correlation functions of the gauge invariant curvature perturbation 
which are observationally accessible in surveys of the cosmic microwave background radiation 
(CMBR) and large scale structure (LSS). C is a popular choice, since it is conserved in single- 
field models once a mode crosses the horizon. However, if more fields are present, any turn 
in the trajectory couples perturbations perpendicular (isocurvature) and parallel (adiabatic) 
to the trajectory [61], which in turn can alter predictions [23]. 

In our setup, we are interested in oscillations of a, which can be the result of a displace- 
ment directly after the tunnelling event or a subsequent event such as a drop or a sharp turn 
in the trajectory (by sharp, we mean 1/Atturn > n^a, in which case damped oscillations with 
frequency nia result [53]). 

There are two effects that can give oscillatory signatures in correlation functions of C- 

1. The direct coupling between perturbations in the fields, which is governed by the rate 
of change in the direction of the trajectory, see e.g. [23, 29, 53]. 

2. The modulation of the background brought forth by cr{t), see e.g. [18-21]. 

The first effect has its origin in multi-field dynamics, that requires a careful investigation 
of perturbations in all fields. To this end, the effective field theory (EFT) of single-field 
inflation in the presence of heavier fields has been developed (not to be confused with the 
EFT of single- field inflation [62-64]), to account for models that are consistent with the 
rational of a low energy effective theory [24-29], see also [22, 23] (some concerns regarding 
the applicability of an EFT raised in [65] were primarily due to different definitions of the 
term EFT, see e.g. the discussion in [22]). The resulting action for perturbations is that of a 
single scalar field with a time dependent speed of sound Cg (t) , which is set by the background 
evolution. While a kinematic basis for fiuctuations [24-29] is appropriate for a slow turn (see 
also [54, 55, 58]), a mass basis is better suited for sudden turns [53]. However, the effective 
field theory tends to break down easily for sudden turns, requiring a true multi- field analysis 
at the perturbed level. Considerable progress has been made in the last years to understand 
signals in the late time adiabatic mode caused by such multi-field dynamics. However, in 
most cases studied analytically, the background evolution is approximated to be smooth, 
i.e. H ~ const is assumed during the turn. 

The second effect takes into account the change in background quantities induced by os- 
cillations of heavy fields, leading to additional oscillatory components to e.g. H{t) or slow-roll 
parameters [18-21]. At the perturbed level, a truncation to the light field is performed, that 
is, any interactions with perturbation in the heavy field(s) are ignored. It is then relatively 
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straightforward to use the in-in formahsm to compute additional features in correlation func- 
tions caused by resonances between mode functions and oscillations of background quantities. 

In a realistic scenario both effects are present. Unfortunately, to our knowledge, the 
added complexity of the problem in presence of an oscillating background and true multi-field 
dynamics for perturbations has so far prevented a comprehensive investigation of observa- 
tional signals (see however recent advances such as the case study in [59]). In the meantime, 
we may focus on either of the two effects to estimate the feasibility of creating detectable 
signals - if it turns out that signals are unobservable for the second effect, one may confi- 
dently constrain oneself to the first one (and vice versa). However, if both effects lead to 
observably large signals, one should follow up with a study combining them, before comparing 
predictions with data. 

In this paper we focus on the second effect, that is we truncate perturbations to the 
inflaton, while ignoring perturbations in the tunnelling field. To apply this logic, any coupling 
of perturbations in the heavy field, Sa, to the fluctuations in the light one, Sep, have to be 
much smaller than existing terms in the equations of motion for 6(p - otherwise the treatment 
is inconsistent. The absence of concrete viability criteria for the truncation in the work by 
Chen et.al. [18-21] led to some lively discussions. Hence, we would like to elaborate briefly 
on the conditions under which such a truncation is viable. 

3.1 Truncating Perturbations to the Light Sector 

Consider the gauge invariant field perturbation 

Q,^S^l + ^tP, (3.1) 

where / G {(j), a}, 4>^ = (/>o- = cr, and ^ denotes the ii perturbation of the metric [66]. The 
equations of motion for Qi are of the form [61] 

Qi + 2,HQi + %Qi + y^MijQj = ^, (3.2) 

J 

where the mass matrix entails the square of the fields bare masses as well as simple combi- 
nations of the background solutions. Diagonalization of the potential's mass matrix Wij = 
d'^W / d(j)id(j)i yields the mass basis of the fields in the potential W[(j), a); for us, this basis is 
given directly by a and (p, since we consider a separable potential. To be concrete, we take 
a quadratic potential for the inflaton in this section 

V{cP)^^ml<P^; (3.3) 

a generalization to other potentials is straightforward. 

Directly after tunnelling, perturbations in both fields should be of the same order, 
Qa ~ Q^, and since Qa decays faster subsequently due to rria ^ m^, we can safely take 
Qa i$ Qcft- Hence, to self-consistently ignore Qa, we need M,^,^ » M^o- as a necessary 
condition. To compare terms in the mass matrix, we take w H^^^., assume slow-roll for 
the inflaton and damped oscillations for a with initial amplitude co- Since perturbations 
in the inflaton need to satisfy the COBE bound, we take ~ \/87r x 10~^, while sixty 
e-folds of inflation requires ~ 16 at the onset of inflation. Thus, to estimate the order of 
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magnitude of terms, we take 



l<Al~^ , l<^l«l^<A|, (3.4) 
\a\ ~ ao , \ct\ ~ aorua , \a\ ~ croml . (3.5) 



We further define 



g2 = _p^^yM<i (3 6) 

which approaches one towards the inflationary regime. To be precise, we have Slaf'^) 



1/V + 1 0.1 at the end of the curvature dominated regime, S{aeq) = l/\/2 at equahty 
and S{af) = l/VWTl ~ 0.995 at the be ginning of the inflationary regime. We first note 



that 

.*^||) (3.7) 

to leading order in slow-roll. The contributions to the mixed mass matrix elements are of 
order 

I ; ; H (ha (ha \ 
M<^.~OUa,0a-^,^,^j . (3.8) 

During the curvature dominated phase H/H^ — )■ — 1, while it is slow-roll suppressed during 
inflation, so that the second term on the r.h.s. is at most of the same order as the first one. 
Further, since \(j)\ <^ \H(p\ due to slow-roll, we see that the fourth term is smaller than the 
first one and we only need to compare the first and third term with m^. Writing H ~ (j)m^/6 
and using (3.4) as well as (3.5), we arrive at the conditions 

c^o<— ^, (3.9) 
rria 

ao « -icpj^ ■ (3.10) 

Thus, the heavier the tunnelling field, the lower the allowed amplitude of oscillations in a if we 
want to neglect the effect of 6a on 5(1). The bounds are strongest in the inflationary regime, 
but get relaxed during curvature domination, since (5 <C 1. If V((j)) is more complicated, 
similar bounds result; we will resort to the above quadratic potential whenever the need for 
numbers arises. 



3.1.1 Exciting Heavy Fields 

The heavy field a is usually displaced from the inflationary valley directly after tunnelling 
[17] ^, but it can also get excited subsequently, for example by encountering a feature in the 
potential leading to a turn as in [23]; it could also run into a sharp drop in the a direction, 
which might be more appealing: consider that after tunnelling, a finds itself on a reasonably 

*This initial excitation is thus similar to the initial excitation of heavy fields considered in [67] . 
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flat plateau of the potential, moving slowly due to large Hubble friction; after some time, but 
still during curvature domination, a encounters a sudden drop, leading for instance into the 
quadratic well with rria ^ rrifp that was used above; the dynamics are thus reminiscent to 
the ones of the waterfall-field at the end of hybrid inflation. We expect that the excitation 
event has a strong effect on the perturbations of the involved field(s). However, if only the 
(T-dynamics are affected, e.g. due to a drop, fluctuations in cp are oblivious to the excitation 
event itself, as long as the conditions in (3.9) and (3.10) are satisfied and Qa is not amplified 
above Q^. A turn in the trajectory is more problematic than a drop, since it actively mixes 
up the roles of the fields, requiring a careful multi-field treatment at the perturbed level [23] . 
Hence, we take as a working hypothesis that the excitation event is of the former type, so that 
fluctuations in (p are unaltered by the excitation of a. It would be interesting to investigate 
such a drop numerically, which goes however beyond the scope of the current article. 

There is a further constraint we need to impose though: since we treat pa as a small 
correction to the background energy, we need to demand pa <C Pk during the curvature 
dominated regime and pa <C <P/2 during inflation^. Using quadratic potentials for both 
fields as before, these conditions entail 



during curvature domination and 



ao«'^l (3.11) 



cro«l'^ (3.12) 
3 nia 



during inflation, which are comparable to the constraints in (3.9) and (3.10). We enforce these 
constraints throughout and, keeping the above caveats in mind, treat ao and its excitation 
time as free parameters. 

3.2 The Curvature Perturbation, Correlation Functions and Mode Functions 

After having justified a truncation to perturbations in the inflaton, we would like to intro- 
duce and comment on quantities of interest, particularly correlation functions of the gauge 
invariant curvature perturbation C, 

(n = (0|C(ki,0)...C(k„,0)|0) . (3.13) 

The two point function defines the power-spectrum Vt^ 

(C'> = ^(2vr)^<^(ki + k2), (3.14) 

which, in single-field slow-roll inflation, is 

Pc = ^. (3.15) 



^The kinetic energy, not the potential one, matters here, since the former determines e = —H/H^ which 
should remain largely unaffected by any sudden changes in pcr and its derivatives during the excitation (for 
example, if kinetic energy in (j) is used to excite a, whose potential energy red-shifts as , H increases 
temporarily); such a sudden change can lead to considerable particle production in Q^, with an accompanying 
bump and ringing in the power-spectrum [68] (similar to the signals in [69] caused by a jump in the mass of 
the inflaton), in addition to the signals computed in this paper (see also [23]). 
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The three-point function defines the bi-spectrum 



p2 

(e) = Sih, k2, h)^{27ry5\ki + k2 + ka), (3.16) 
11 

i 

which is commonly characterised by a dimensionless shape function S; we use the conventions 
of [18] for ease of comparison, see App. A. Translational invariance requires that the three 
momentum vectors form a triangle, and rotational invariance implies that only the amplitudes 
of the momenta enter (ki = |kj|). The shape of S refers to the dependence on the shape 
of the triangle formed by the three momenta (i.e. the dependence on ki/k^ and k2/k^ while 
keeping K = ki-\-k2 + k^ fixed) ; the running of S refers the dependence of S on K for a fixed 
shape. 

What determines The equation of motion for adiabatic perturbations can, via field 
redefinitions, always be rewritten as a simple harmonic oscillator with a time dependent mass 
(set by the background evolution), see for example the review [66]. The variable satisfying 
this equation is commonly referred to as the Mukhanov variable v^- 

On small scales [k ^ aH), the mass term becomes negligible compared to the gradient 
term k'^Vk-, and simple oscillatory solutions result for Vk- Hence, the presence of spatial 
curvature is irrelevant for the evolution of on small sales. The amplitudes of the two 
independent oscillatory solutions are usually set (after quantization) by imposing the Bunch- 
Davies vacuum state. 

Not surprisingly, the oscillations in are directly related to oscillations in C,: if we 
define the mode function of Q by 

Uk = j d^xC{t,x)e-'^-'' , (3.17) 

it is easy to show that a Bunch-Davies state for translates to 

Uk « ,j , e-'"^ . (3.18) 
a-y/4ebgrfc 

These oscillations on small scales are of particular interest to us, since they can resonate 
with oscillations in background quantities. In the above, we expanded the full solution in 
the limit k 3> aH, which is a good approximation for all cases discussed in this article. Once 
scales cross the horizon, one can show that ^ is frozen; subsequently, ^ sets the amplitude of 
temperature fluctuations in the CMBR, enabling a direct probe of early universe physics. 

Before we compute correlation functions of we need to address a crucial point: the 
initial state of fluctuations after the tunnelling event in open inflation is not the Bunch Davies 
state; a detailed discussion of the state for single and two-fleld one-bubble open inflation can 
be found in [9]. 

The deviation from the standard Bunch-Davies state arises due to the mass difference 
outside the bubble wall and is parametrised by 1 — y with 

T(2-ip) e-'^PBr, 

y = -^ -T , ^_ +C.C.. 3.19 

r(2 + ip) 2 cosh vrp ap 

Here otp and /3p are coefficients relating the mode function of the true vacuum inside the 
bubble to the mode functions of the false vacuum outside the bubble, and p labels momenta 
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(for more details, see [9]). The key observation to be taken away is that on scales smaller 
than the curvature scale (sub-horizon modes with large p), Y tends to 0, and the state can 
be approximated by the Bunch-Davies one. This agrees with our intuition (and the equiv- 
alence principle) that given a small enough patch, spacetime inside the bubble approaches 
Minkowski. 

Hence, we proceed with the Bunch-Davies state as long as we are concerned only with 
modes inside the horizon. We would like to emphasize that any (small) deviation from the 
Bunch-Davies state inside the Horizon might slightly affect the amplitude and frequency of 
Uk, but not the presence of high frequency oscillations that are key to our upcoming analysis. 

3.3 Overview of the In-In formalism 

The in-in formalism was initially introduced by Schwinger [30] and used more recently in a 
cosmological setting by Maldacena [70] and Weinberg [71] among others (for recent reviews see 
[31] and [32]). It can be used to compute correlation functions of the curvature perturbation 
C, for example the power-spectrum and the shape function of the bi-spectrum. Here, we 
provide a brief summary following [18], primarily to introduce our notation. 

Starting from a Lagrangian, one may perturb it around a homogeneous solution of 
the classical equations of motion; quadratic and higher order terms provide the dynamics 
of perturbations. Conjugate momenta are defined from the interaction Lagrangian, enabling 
the computation of the Hamiltonian in the perturbations. By gathering all terms higher than 
2'nd order, one arrives at the interaction Hamiltonian Hj. Solutions of the classical equations 
of motion for the quadratic part of the Hamiltonian in Fourier space are commonly called 
mode functions, which are normalized according to a Wronskian (quantization) condition 
with initial conditions corresponding (usually) to the Bunch-Davies vacuum. A correlation 
function is then computed as the nested integral expression 



iCit)) 







Texp 



H^dt 



C{t) 



Texp 



H^dt 



to 







(3.20) 



where T is the time-ordering operator (T the reverse one). After expanding H to the desired 
order, one can perform contractions, sum over them and perform a time-ordered integration. 

Our aim is to compute the two and three-point correlation function of C- The Hamil- 
tonian H (not to be confused with the Hubble parameter, which does not play a role in 
this section) is a functional of the field (p and its conjugate momentum it. It determines the 
system's evolution via 



i[HicPit),7Tit)),<Pi^,t)] 



^(x, t) = i [H {(Pit), nit)), 7r{K, t)] . (3.21) 



We expand the field in terms of a time-dependent background (denoted by a bar) and per- 
turbations, i.e. (pi:x.,t) = (pit) + 6(pix,t). The background evolves according to the classical 
equations of motion, 

- fyi-l 

m = m = (3.22) 

where 7i is the Hamiltonian density. Since the background commutes with everything (it is 
just a complex number), the canonical commutation relations are transferred to the pertur- 
bations alone, 



[6(Pi^,t),57riy,t)]=iS'i^-y), 
[6(Pi^,t),5cPiy,t)] = [57ri^,t),67Tiy,t)] = 0. 



(3.23) 
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Taylor expanding the Hamiltonian around the background solution yields 

/BT-I I' BT-l 

+Ho {6^{t),6TT{t),t) + Hi {5<P{t), 6Tr{t),t) . (3.24) 

Here, Hq contains quadratic terms in the field, which determine the leading evolution of the 
perturbations, and Hj contains all higher orders. Combining these results, we see that the 
equations of motion for the perturbations are determined by Hq+Hj. The point of separating 
the non-linear Hamiltonian in two components is to determine first the kinematic evolution 
as driven by Hq, and then perturbatively compute the effects of interactions induced by Hj. 
For a generic operator Q in a state the expectation value of Q is [18] 

{n\ Q (J0(x, t), <57r(x, t)) \n) = {n\ F-\t, to)Q {6^^ix, t),5TT\x, t)) Fit, to) \n) (3.25) 

with to in the far past and 

F{t,to) =rexp J Hi{t)dtj . (3.26) 

Hence, for Q = C (3.25) reduces to (3.20). 

Note that after quantisation, the introduction of mode functions, and application of 
T and T, the creation and annihilation operators kill the respective left and right vacuum 
states. Hence, the remaining terms are those where the fields are contracted away. Further, 
one can use Feynman diagrams to keep track of the necessary contractions, while keeping 
in mind that the proper vacuum state has to be prescribed in order to compute an explicit 
integral. 

3.4 The Power- and Bi-spectrum 

To leading order in Ae = eosd and A(e?i) = {eT])oscij the interaction Hamiltonian can be 
expanded to [19]^ 

-Ho = ahoC^ - aeoidCf , (3.27) 
ni ^ -a^AeC^ + aAe{dCf , (3.28) 

ni « -^a'Aiev)eC + -.. ■ (3.29) 

We only kept the most important term in (the one with the most time derivatives acting on 
slow-roll parameters, as in [72, 73]). Thus, the leading correction to the two-point correlation 
function is given by 

A (C'> = (ON j adrd^x [Ul CkiCfe] |0) . (3.30) 

Using the definition of the power-spectrum in (3.14) and the unperturbed power-spectrum 
of a single-field slow-roll model in (3.15), one arrives at 

APt /"Tmax / 9 \ 

-^(A:) = 2i dra^Ae (u\ - k^ulj + c.c. . (3.31) 

C '^min 



These results were derived in a spatially flat universe, but they remain valid in the regime of interest if 
curvature is present. 
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Similarly, the leading order contribution to the 3-point correlation function is 

''"max 

<C'> = {0\i I adrSx ["H^ ('(kikskg)] |0) 

du 







dr 

+2 permutations + c.c. , (3.32) 
so that the shape function in (3.16) becomes 

S (fei,fc2,fc3) = I I dra A{er])Uf,^Uk^-j^ + 2 perm. + c.c. . (3.33) 

Analytic approximations for (3.31) and (3.33) in inflationary regimes are given in [19], which 
we recover (and slightly extend) in section 4.2. The corresponding results in the presence of 
curvature are provided in Sec. 4.1 and Sec. 4.3. 

Substituting the mode function into (3.31) and (3.33), and incorporating the complex 
conjugates, the expressions simplify to 

^'PC^,^ r-"^" //(t)2-2A:2 , , A Ae 



and 



{k)= sin(2A;r) + 2/(r) cos(2A;r) —dr (3.34) 



S {kiMM) = kik2k. I ^ r™^|^(_2i^cos(ifT)+6/(r)sin(i^r)) dr , 

(3.35) 

with /(t) from (2.22), encoding the time-dependent background evolution of the amplitude 
of Ufc, and we defined rescaled wave- numbers 

r ki ~ k ~ K kT + ko + k-i , , 

"'eq "-eq "-eq "'eq 

where 



(3.37) 

is the wave-number corresponding to the Hubble horizon at equality^, pk = p^', for ease of 
notation, we dropped the tilde 

ki ^ ki , etc. (3.38) 

In the remainder of this article, all wave-numbers are normalized in this manner. 

The terms containing /(r) in (3.35) and (3.34) are missing in [19], since they are sub- 
leading during inflation, but they play an important role if curvature is important. The 

'^In all plots, we choose appropriate approximations for H in the three regimes, e.g. during the curvature 
dominated regime we use H ^ \/pk/3- 



keq — Oj(teq) H (^teq) 
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pre-factor fciA;2^3 in (3.35) is a due to the chosen normahzation, see Appendix A; thus, we 
plot S/{kik2ks) in figures. 

The integration hmits above must be adjusted to the boundaries of the region under 
consideration - this is crucial for the transmission of a signal, see Sec. 4.1, 4.2 and 4.3. As 
discussed in section 3.2, we take the initial state to be the Bunch-Davies vacuum. In order to 
compute the integrals (3.34) and (3.35), we note the high-frequency oscillations of the mode 
function (3.18) on sub-horizon scales and the oscillations of the slow-roll parameters induced 
by the massive u-field, whose frequency increases with the mass of the field. They allow the 
leading order behaviour of the integrals to be analytically approximated via the stationary 
phase approximation. For details on this asymptotic method, see e.g. [74], chapter 4. The 
relevant formula is 

/oo 
i;{z)e'^^'Uz ~ V(^*)e'(^("*)+''^''(^"("*))^V2^/I^"(^*)l, (3-39) 
-oo 

for F{z) — )• oo and ip{z) sufficiently well-behaved (smooth and without high-frequency os- 
cillations) as well as F'^z^,) = 0. Away from this equilibrium point, the integrand oscillates 
rapidly so that the contribution to the integral cancels out approximately over one period. 
If the equilibrium point falls within the limits of integration in one of our integrals, a signal 
is generated; if it falls outside, the contribution vanishes. This asymptotic approximation 
introduces a step function at the boundary, i.e. when the equilibrium point coincides with a 
limit of integration, so that no remnants of a signal can seep through it; in reality, a con- 
tinuous transition occurs. Since we rely on this analytic approximation, we expect to obtain 
results with sharp boundaries and possible discontinuities. 

We focus on analytical approximations in order to gain insight into the resonance mech- 
anism. Furthermore, due to the highly oscillatory integrand, a numerical computation is 
challenging and usually of limited practical use, since no direct detection of oscillations in 
the power or bi-spectrum has been made yet. 



4 Oscillatory Signals in Correlation Functions 

We have now the tools at our disposal to compute the signals in the power-spectrum (3.34) 
and bi-spectrum (3.35) brought forth by the oscillating tunneling field a. We discuss the three 
regimes introduced in Sec. 2.2 consecutively: curvature domination in Sec. 4.1, inflation in 
Sec. 4.2 and the intermediate regime in Sec. 4.3. 

4.1 Curvature Domination 

In a curvature dominated universe we can ignore the inflaton field, and solve the Friedmann 
equation (2.4) to 

at , . 

We start our discussion of the curvature dominated regime at a(to) = o,'^^ = 1 with to = I 
so that a = t. Further, for computational convenience, we work with conformal time r, 
normalised to the time of equality, 

^0 ,„ / ^eq 
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In other words, we convert the entire region of interest (plus a neighbourhood) to the negative 
real axis via the map 1 1— r : (0, teq) ^ (—00, 0). This leads to the following relations for the 
conformal time^, the scale factor, the background Hubble parameter, and the background 
slow-roll parameter 

T = In , a(r) = y^p^e^ , i^bgr(i) = ^ , ebgr = 1- (4.3) 

The solution of (2.5) for a is given in terms of Bessel functions with the asymptotic limit, 

f^cTof-^) I sin (mat + a) + cos (mat + a) | , (4.4) 

\toJ V 8m<^t J 

for mat ^> 1 and a a phase set be the details of the initial excitation mechanism; the latter 
is of little interest to us (we often set a = for simplicity). The subscript "0" indicates the 
moment of the initial excitation, e.g. the tunnelling event. The maximal oscillation amplitude 
is constrained by (3.9) and (3.10) (for larger amplitudes, multi- field effects become important 
at the perturbed level, which can not be recovered in our approach). The resulting oscillatory 
contribution to the Hubble parameter becomes 

i?osci = -^(^) 'sin(2(m.t + a)) , (4.5) 
and the leading oscillatory components of the slow-roll parameters read 

eosci = T^ - cos(2(m,t + a)), (4.6) 
4-Hbgr V*o/ 

(e^)osci = (^) ' COS (2 {mat + a)) . (4.7) 

Using the above in the power-spectrum correction (3.34) and the shape function (3.35) leads 
to 

^''^w= .,i'"K r"(-v' (4.8) 



n ifl^ Co) 7,,,.,. V'o 

cos (2 {mat + a)) ( -j^ — sin(2A;r) + 2 cos(2A;r) j dr, 

2 4 

S{kiMM) = fcifc2fc3AA (4.9) 

16-f^bgr(io) 

{2K cos{Kt) - 6 sin(Kr)) cos(2(m^t + a))dT, 

■ mill 

where we recall that all ki are rescaled with A;cq, we used K = (/^i + /c2 + ^3) and defined 
®Thus, conformal time is identical to the number of e-folds normalised to acq, r = A'' = In (a/ acq). 
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Since ai^bgr = 1 during the curvature dominated regime we have M = 1. After substituting 
the integration boundaries of the curvature dominated regime, 



r ■ - t'-''^ t - T^^'^ (A 11) 

'mm — 'j ) 'max — 'j V^'-"-^/ 



from Sec. 2.2, and applying the stationary phase approximation, the power-spectrum becomes 
for high frequencies of the arguments of the trigonometric functions (large k and/or ttt-o-) 



+ COS (2fc(l -a;(A:)) + 2a + ^) y (4.12) 
where the frequency is set by the 

.;(fc)^lnf '^ ) . (4.13) 
Similarly, the shape function can be approximated by 



S (ki, k2, k3) « kik2k3M f COS {k (1 - u^iK/2))) + 2a + ^ 

^ sinfif(l-6^(i^/2))) + 2a + ^)y (4.14) 



K \ ' ^ ' 4 ^ 

We plot the resulting oscillatory corrections to in Fig. 1 for exemplary choices of the 
excitation time t\ (we allow t\ to differ from to? in which case the excitation is not a result 
of the tunnelling event), initial amplitude ctq and mass mo-, making sure that gq satisfies 
(3.9) and (3.10). Evidently, the oscillation amplitude is proportional to cJq with a decaying 
envelope at large wave numbers, caused by the /-terms in (3.34), with /(r) defined in (2.22). 
Corresponding oscillations are present in S, plotted in Fig. 2. However, its envelope has 
a minimum due to the decaying sine mode induced by the fir) term. For comparison, the 
corresponding cosine mode in the power-spectrum has little effect, because it is sub-dominant 
over the entire k-domain. 

For both, power-spectrum and bi-spectrum, the cut-off in the signal is determined by 
the length of the curvature dominated period between 

l = o!i^<a,<a< af^ = ^ . (4.15) 

Since the resonating wave- number leading to the oscillations is given by fc = arua, the 
boundaries of the scale factor translate to the interval 

m^<k<—^^, (4.16) 

(k) 

for our rescaled wave-number if Oi = Oq (and a shorter interval if the oscillations start at 
ti > to), where k stands for K/2 in the shape function case. 

It is instructive to compare the period of the cr-field oscillations T = 27r/m„ with the 
available time tf-to = {A'^Pcp/^y'^^'^ - 1- If a fj-oscillation can't be completed before the end 
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Figure 1: The power-spectrum in the curvature dominated regime for different masses of the cr-ficld 
plotted over the physical wave-number fc, normalised by the horizon scale at the time of equality 
kcq = a(icq)^^bgr(icq), with a = 0. Tlic time of excitation ti is taken to be after the tunnelling event. 

of the regime, a resonance with the vacuum oscillations is impossible and no power-spectrum 
signal is generated. However, since the individual resonating modes ki involved in the bi- 
spectrum calculation are necessarily smaller than the corresponding power-spectrum modes, 
they have a chance to resonate with the u-field, leading to a signal in the shape function, 
even if none is present in the power-spectrum. Given a few oscillations, a small resonant 
running arises. 

4.1.1 A Case Study and Observational Prospects 

It is instructive to consider a concrete inflationary setup, particularly to estimate the maxi- 
mally attainable signals in the power- and bi-spectrum. 

As in Sec. 3.1, we use a simple inflationary model with V^(0) = m?^(fp' ~ -v/Svr x 
10~^ and w 16 in the slow-roll limit, so that around sixty e-folds of inflation result after 
curvature becomes sub-dominant^ and the normalization of the power-spectrum is satisfied. 

®That is, we put the observational window of the CMB around fcoq. Since signals are generated deep inside 
the horizon during curvature domination, relevant modes can cross the horizon during inflation and he in the 
observational window. 
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(c) m^r = 15v/p0/3, CTo = 0.1 

Figure 2: The shape function in the curvature dominated regime for different masses of the cr-field 
pfotted over the physical wave- number k, normahsed by the horizon scale at the time of equaUty /coq, 
with a = 0. The time of excitation t{ — (1500-\/ p^/i)^^ is taken to be after the tunneUing event but 
well before the end of the curvature dominated regime. 



The curvature dominated regime ends at 

3.0 X 10^ (4.17) 



where we used A = 10, which implies an interval in (4.16) of 

777 

m^<k< —0.015 (4.18) 
where a signal could be generated, and an oscillation frequency in (4.13) of 

For which ma can a signal be generated in the power-spectrum? Allowing for at least 
one full oscillation of a in the curvature dominated regime leads to the condition 2ir/ma < 
t{ — to = af^^ — a[^^ , which becomes 

^ > « 2.A^ ^ 410 . (4.20) 
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Thus, the tunnelhng fields mass needs to be bigger than 

^ m,, (4.21) 

considerably heavier than the mass of the inflaton. This is not surprising, since otherwise a 
is frozen during the curvature dominated regime due to Hubble friction, just like the inflaton. 

Assuming rn-o- > m™'"^, how big of a signal can be generated? Since 5 < 0.1, see 
eqn. (3.6), and <C rria, the strongest constraint on the oscillation amplitude of a stems 
from (3.9), that is 

<r«"-^\. (4.22) 

which needs to be satisfied throughout the oscillations of a. Violating this bound does not 
mean that signals are absent, but that they can not be computed in the truncated single- field 
description at the perturbed level, which we employ. 

Since 5{t) oc a{t) and the oscillation amplitude of o" redshifts as a~^/^, the bound remains 
satisfied, if it is satisfied initially. Hence, we need to impose 

c7o«^^, (4.23) 

at the time of the initial excitation. The farther the excitation lies in the past, the smaller 
is the observable signal. To provide a conservative upper bound on amplitudes, we consider 
an initial excitation of a at exactly one oscillation period before the end of the curvature 
dominated regime, 

T= — =t{-ti = ^ - Oi , (4.24) 



which determines ai; thus, we arrive at 

5i = ) , (4.25) 

with (^f'^ = + 1 ^0.1. Such a single period in a generates a signal close to the upper 

boundary of the interval in (4.18) only (i.e. in Fig. 1 before the upper cut-off on k). 

Using 3i?^gj.(to) ~ ^{4')/Si ^ ^(0) ignoring factors of order one, we can put an 
upper bound on the amplitude of the power-spectrum correction in (4.12) 

APr I , , 

< 5i — (4.26) 



Pc_ (j)^ 



1 f ^^27) 



m 



0.01 ( - 0.9997 ) , (4.28) 



and the bi-spectrum in (4.14) 



S ...oml 1 ..fAP^ ^ 



C 



«-^'^Si4=-^(^) (4-29) 
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where M = 1. Thus, a correction at the percent level in the power-spectrum is possible if 
the initial excitation happens close to tf and rricr is large enough; of course, any mass should 
be smaller than the Planck mass, 

ma<.l, (4.30) 

yielding AP^/P^ ^ 5 in our toy model. In the plots of Fig. 1, we chose earlier times to 
cover several oscillations, and masses well below the Planck mass, resulting in weak signals. 
(To satisfies the bound in (3.9) in all plots; in these cases the resulting amplitudes are much 
smaller than the above bound. However, even if the correction to the power-spectrum is 
about 10% (possible to generate with close to the Planck mass and tuning of ti, while still 
marginally allowed by current data), the bi-spectrum remains dauntingly small, S/ {kik2k^) ~ 
0(1). Further, as mo- is increased, the frequency uj in (4.13) grows too, complicating a 
detection in a noisy data set. 

Let us estimate the amplitude of the power-spectrum correction if the initial excitation 
is due to a displacement of a directly after the tunnelling event. Using the same line of 
thought as above with ai = oq = 1 leads to 

« (4-3.) 



Pr 



m, 



« 3.3 X 10"'' — ^ . (4.32) 

Unfortunately, such a signal is currently unobservable: even if rua = 1, we have AP^/P^ ^ 
0.0016. 

Lastly, let us have a brief look at current observations. The amplitude of any oscillations 
on top of the power-spectrum of the WMAP satellite has to be lower than about 10% [21, 75- 
77], depending on the frequency. Such oscillations in the power-spectrum are conceivable, 
but do not improve the theoretical fit significantly [76, 77]. Using the flat sky approximation 
cLa = l/k, where dA = 13.7 Gpc, and the pivotal scale /spivot = 0.027 Mpc^""^, we can estimate 
the maximal frequency resolution to 

, , 1162 

w < , 4.33 

s 

where 

s = ^ , (4.34) 

and we required that an oscillation stretches over at least two /-modes in order to be observ- 
able. Realistically, some binning of the data is necessary to reduce the noise. Combining 
10 Z-modes into a single bin, the maximal frequency is reduced by a factor of 10; hence, we 
define 

116 

= — • (4.35) 
We can tune s, and thus the position of this window 

feobs e (0.0054, 3.2) s, (4.36) 
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by moving the initial value of the inflaton and/or its potential after the tunnelling event, as 
long as the subsequent inflationary phase last long enough to dilute curvature below current 
observational limits. Naturally, we are primarily interested in s ~ 0{1), since otherwise any 
trace of the curvature dominated regime is pushed out of the observable range. 

Fortunately, the frequency in (4.19) is only logarithmically dependent on k and nifp/mij: 
taking s = 1, k = 0.0054 and rria = m™™ yields the estimate 

1^1 < 9.4, (4.37) 

well below u^^. 

To summarize, the main challenge to observe an oscillatory signal in the power-spectrum 
generated during the curvature dominated regime is the generically small amplitude: an 
observation would indicate an excitation shortly before teq with a near maximally allowed 
amplitude and a large m„. However, the frequency of oscillations is commonly low enough 
to pose no additional challenge. The generated bi-spectrum is proportional to AP^ / and 
thus unobservable with current experiments. 

4.2 Inflation 

The inflationary regime has been discussed in [19], whose results we recover and extend 
slightly. Neglecting the contribution from curvature, we consider a simple power-law infla- 
tionary regime, 

acxtP. (4.38) 

Conformal time is given by 

It , , 

T=- 4.39 

{l-p)a 

As p — )• 1, a conformal singularity is approached, indicating the need to treat the transition 
region to curvature domination separate case. 

Using the Hubble parameter -ffbgr = p/t, the solution of the cr-field's equation of motion 

is 

fj (To I — I sin (mfjt + a) -\ cos (niat + a) . (4.40) 

\toJ V J 

The initial amplitude may be set during inflation or in a prior regime, in which case to 
coincides with the beginning of the inflationary era. We only focus on the latter case in this 
article, keeping in mind that a generalization to the former one is straightforward. Thus, the 
leading order oscillatory contribution to the Hubble parameter is 



Hosci = -^Y^\y-j sm{2{m^t + a)), (4.41) 
so that the leading order oscillations of the slow-roll parameters become 

-3p 

= 777?^ ( f ) cos (2 {m^t + a)) , (4.42) 





to) 






2 4 




-"bgr 


\to 



'3p 

(e^)osci (t) = ( f ) COS (2 (m^t + a)) . (4.43) 
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Inserting the above in the power-spectrum correction (3.34) and the shape function (3.35) 
leads to 



(k) 



alml 



4^bgr (*o) ebgr 

2 



-3p+2 



cos (2 {rrifjt + a)) 



{AAA) 



k{l — pYt^ 



2k sin(2/cT) + 



2p 



{1-p)t 



cos(2A;r) I dr , 



and 



2 4 



S{ki,k2,k3)=Mkik2k3 ^ 



3-3p 



(4.45) 



X cos(i^r) 



3p 



[1-p)t 



sin(il'r) I cos(2(mo-t + a))dT . 



respectively, where we used the definition of the normalization constant M from (4.10), 
reading 



AA(ieq, to) 



to 



3(p-l) 



(4.46) 



for power-law inflation. Based on the stationary phase approximation, the asymptotic be- 
haviour of the power-spectrum correction 



Pr 



(k) 



¥ 1 



^ f2k\ 



4 ebgr ° \H^gr{to)J \krj 



(4.47) 



bgr 

2m'i 



+ 



kv J 

1 

1 f2k\ p -ffbgr(io 



) ■ )sin(^(A:)^ + 2a+^ 



. 2k TT 

cos ( ijj[k)- — h 2a + — 



and the shape function 



5(/ci, /c2, fcs) ^ kik2k3M 



TT 1 



8 Ebgr ° V-f^bgr(io)/ V^r 



5 



2p 



x|c„s|.(A72)- + 2a.+ -)+^(- 



(4.48) 



sin ( oj {K/2) ^ + 2a + ^ 



is again oscillatory, as in the curvature dominated phase, albeit with a different frequency 



vo{k) = p 



P 



1 — p 



2m„ /2k\ 



Hhgr{to) \kr J 



(4.49) 



and pre- factor. Following [20], we denoted the first mode that starts to resonate as soon as 
the cr-field is excited by 



kr = 2maao ■ 



(4.50) 
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Figure 3: Power-spectrum corrections during the inflationary regime with Cbgr = P^^ ■, to — toq, 



a — 0, and kr — 2m„aQ for different values of nic / Pify/i and p 



We note two minor differences as compared to [20]: firstly, we retained the /-term 
contributions to the power-spectrmn corrections, providing the two sub-dominant terms in 
(4.47). Secondly, just as in the curvature dominated regime, we retain a dependence on 
^1^2 ^3 of S stemming from the first term in (4.48) tracing back to our chosen definition of 
the shape function. Retaining the dominant term only, we fully recover Chen's results. 

The expressions (4.47) and (4.48) are valid during inflation 



di < a < 



with af from (2.19) as well as B = 0.1, and 



which translates to the interval 



< k < 



m. 



(4.51) 



(4.52) 



(4.53) 



where a signal could be generated. Here k stands again for i^/2 in the shape function case. 
The boundary is only important at the low end of the spectrum, since the exponential factor 
sends the upper boundary to unobservably high wave numbers. Physically, the cr-oscillations 
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(c) m<,/yp^ = 100, p = 20 

Figure 4: Shape function during the inflationary regime with, Cbgr — P^^, to ~ tcq, cto = 0.0001, 
a = 0, and kr = 2mo-ao for different values of m„/ and p. 



resonate only with modes that haven't crossed the horizon yet, so that we only need to be 
careful about the lower bound of the interval in contrast to the curvature dominated regime. 
If the excitation takes place late during inflation, the signal is shifted out of the observational 
window. For a detection, the amplitude of the oscillations is of prime interest. We saw that 
the amplitude of oscillations caused by an initial displacement of a after the tunnelling 
event is too small at the end of the curvature dominated regime to be detectable in current 
experiments (within our computational framework). Hence, if an oscillation is detected, it 
must have its origin in a subsequent excitation event. Further, the scale dependence of the 
signal's amplitude is model dependent [20]: changing the slope or the width of the slow-roll 
valley can give qualitatively different envelopes of the oscillatory signal. 

For the quadratic potential of cr, the amplitude of the bi-spectrum envelope decreases 
during inflation (Fig. 4), whereas it increases during the curvature dominated regime (Fig. 2) 
for all sets of parameters. Therefore, we expect to find a local maximum during the inter- 
mediate regime, which could act as a tell-tale sign of open inflation, if it were observed. 

The frequency uj{k) of the signal is again set by the inverse of the scale factor. Thus, 
if the amplitude is observable, the frequency of the oscillations acts as a direct probe of the 
scale factor's time dependence, as can be seen in Fig. 3 and 4. We refer the interested reader 
to [19, 21, 75-77], where signals generated during inflation and observational prospects are 
discussed in more detail. 
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Before investigating the intermediate regime, we would like to comment briefly on the 
allowed amplitudes for ao during inflation. Since (5 ~ 1, the constraints on ao in (3.9) and 
(3.10) read 

ao < ^ , (4.54) 

2 

ao<^</., (4.55) 

where we used again V^ip) = m^cf^ /2 as a prototypical inflationary potential. Some of the 
concrete case studies in Chen's work come dangerously close to these limits, particularly if 
rua is large. For example, a case study in Sec. 6.3 of [19] (see also [21]) yielding a 10% 
correction to the power-spectrum requires (leaving out all dimensionless factors of order one) 
(3 ~ (/)^/(m^(TQ) ~ 10^^ during slow-roll, e ~ (jP'/H'^ ~ 10^^, which corresponds to an initial 
oscillation amplitude of do ~ W'^H/nia ~ 10~'^ (pm^ / rrirj , where we used again our toy model 
with ~ 16 in the last step to enable a comparison with our limits. While the first bound 
in (4.54) is marginally satisfied (one has to keep in mind that all dimensionless factors of 
order unity are ignored, which could easily lead to a factor of 10), the second bound in (4.55) 
becomes stronger for rriu/m^ > (/) ~ 16 and is therefore violated easily for heavier masses. 

We conclude that if a detection of oscillations at the few percent level in the power- 
spectrum were to be made, and one wishes to attribute this signal to the oscillations of a 
heavy field during inflation, one should go beyond the current computational scheme and 
include perturbations in a throughout. 

4.3 Intermediate Regime: Comparable Contributions 

In the intermediate regime, the contributions to the energy density from spatial curvature 
and the infiaton field are comparable. Since this interval is brief, Hubble friction is high, 
and a evolves much faster than (p since rriu ^ rricf,, we model the infiaton contribution as a 
constant in this section, 

Ptf, const. , (4.56) 
so that the solution of the Friedmann equation (2.4) is 

a{t) = ^^^^(yg3^) ^^^^^^ 

and the expression for conformal time r takes the form 

r = ln(^tanh(^^^^^. (4.58) 

Note that the units of time are normalised differently - with respect to the moment of equality 
rather than tunnelling. 



' eq 



-arcsinh(l) , (4.59) 



in order to map t G (0, oo) i— )• r G (— oo,0). Using (4.57), the background Hubble parameter 
becomes 

^bgr = y^coth(y^t) . (4.60) 
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Given this background, the equation of motion for a yields no simple analytic solution. In 
order to proceed, we expand the Hubble parameter around the time of equality 



"eq 



arcsinh(l) , 



so that 



■gr 




\f2 + arcsinh(l) 




+ 0((t-teq)2). 



(4.61) 



(4.62) 



Retaining the linear term is sufficient for our purposes, considering that i/bgr acts merely as 
a friction term and we perform an asymptotic expansion later on. 

As before, in order to solve the equation of motion for the c-field, we need to specify 
initial conditions. Technically, they are prescribed by the CdL instanton solution of the 
tunnelling event, dynamically evolved to the neighbourhood of fgq with the solution of the 
equation of motion during the curvature dominated era given by eqn.(4.4). However, provided 
that the frequency of oscillations in the curvature dominated era is large enough (which is 
the case), we can always choose to begin our considerations at the peak of on oscillations at 
to, where a = and = 0. Further, we learned that an initial excitation directly after the 
tunnelling event is too weak to be observed; hence, we envision a subsequent excitation in the 
curvature dominated regime and treat gq as a free parameter, only subject to the constraints 
in (3.9) and (3.10). 

Given these initial conditions, g can be expressed in terms of Hermite % and confluent 
hypergeometric (Kummer) functions iFi as 



G{t) = Go 



\/3 4 + \/2Jp0/3tcq iFi 1 



6/90/3' 2 

y|(t-teq)y^W3-\/3] - An ^2 



^^(2^/2 + y^teq)') 

^/3, 



3P0/3 



(4 + J2p0/3t 



cqy 



^/?,(^+^2p^/^U^ iFi [l 

( Vs 



3 3 



6p0/3'2'8 



(2^/2 + . W3tcq)' 



P0/3 

iFi 



^ (4 + ^2/50/3teq) 



V3, 



(4 + J2p^/^t 



cqy 



1 3 



^;-;-(2^/2 + .W3t, 



6/9^/3' 2' 8 



(4.63) 



Since we are primarily interested in the leading order oscillation frequency and the am- 
plitude's evolution, let us approach the problem from a more pragmatic perspective. The 
equation of motion (2.5) is that of a damped simple harmonic oscillator with mass mo-, 
justifying the ansatz 

g^gq(^ cos{m^{t - ti)) , (4.64) 

which solves (2.5) to order O {mcr/tcq)- Comparing to the full analytic solution, we find that 
the best fit for the exponent is p ~ —1.5. This exponent depends only weakly on rria, but 
changes somewhat with teq (i-e. it depends on y^p^/S via = P(p/3) and the size of the 
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Figure 5: Power-spectrum correction in the intermediate regime for different masses of the (T-field, 
with (To = 0.01, and A — 10, where k^q is the horizon scale at the time of equahty. The power law 
index p is determined by TOg. ■ 



neighbourhood we are considering, ti — t[. Given the notation of Sec. 2.2, the initial and final 
time are 

/ 3 / 3 

ti=J — arcsinh(l/^) , t{=J — arcsinh(l/i?) . (4.65) 
\l P<l> \l P4> 

Since we chose ^ = l/i3 = 10, we have 

/ 3 

if-ti«2.9 J— ss3.3teq- (4.66) 
V P't> 

From here on the computation proceeds as usual. The oscillatory component of the 
Hubble parameter is 

^osci = t(t] sm(2m.(t - 1,)) , (4.67) 

and for a large enough ratio mo-/ y^p^/S > 15, the dominant oscillatory component of the 
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slow-roll parameters take the form 

_ al p ( t 



cos {2m„{t - ti)) , (4.68) 



N)osci = ~^t(t] cos {2m At -U)) , (4.69) 



leading to the power-spectrum correction 

APc ^ «P flV'" ^^^^^^W^ eos(2n..(t-tO) 



Pc Qp^/^UJt^^^ \tj cosh(v^t) 

.2/^^ _ 01,2 



/ la ^ sin(2A;r) - 2tanh(r) cos(2A;r) j dt (4.70) 



and shape function 

Qd. I. Kfl-l-l. 4<V^rP ( t Y^-' Sinh3(y^t) 

S {ki,k2,k3) = J\fkik2k3 i i \ 



24i7bgr(io)tr\/P073 -^--in cosh(7V3t) 
2K cos{Kt) + 6 sech I J sin(Kr) j cos (2m^(t - ii)) . (4.71) 



with in (4.10), which becomes 



2^2 x3/2 



where we used our approximate background solution in (4.57) and (4.60) as well as A = 10 ^ 
1 from Sec. 2.2 in the last step. The asymptotic behaviour of the power-spectrum is given by 



Pt; ^ ' ^ 12arcsinh(^-i) ° \ ml J 



/arcsinhf^^^\ 



2p-l 



arcsinh(yl ^] 



2^ (^1 + - l) sin (^2kojik) - 2m^ti - ^ 



+2k ^1 + ^ ^J^ ^ ' cos (2kuj{k) - 2mati + j > (4-73) 

with frequency 



uj(k) = ( , arcsinh [ VP4>/ _\ _|_ g^j-^g^j^]^ | | ] . (4.74) 
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Figure 6: Shape function in the intermediate regime for different mass of the (j-field, with (Tq = 0.01, 
^0 — tcq so that a (to) ^^bgr(io) = where k^q is the horizon scale at the time of equahty. The power 
law index p is determined by TOo- . 



Similarly, the shape function becomes 



S{ki,k2,k3) 7^ Mkik2k^ 
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(4.75) 



These signals, plotted in Fig. 5 and 6 for a few exemplary parameters, show a mixture 
of familiar features: firstly, the signal is bounded in the k-domain on both sides, 



(4.76) 



set by the duration of the intermediate regime (we chose A = 1/B = 10, and k stands for K/2 
in the shape function case) . Note that signals from the curvature and inflationary regime are 
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continuously connected in momentum space by the intermediate interval, with later times 
corresponding to higher wave numbers. 

If tf — t[ < T = 27r/m(j, no signal is produced in the power-spectrum (a resonance 
requires at least one full oscillation of a), which translates into a lower bound on given 
the duration in (4.66). 

Two unexpected features are evident in our analytic results: 

1. The amplitude of the oscillations in the power-spectrum exhibits a local minimum 
and increases above a threshold wave number; we expect the amplitude to decrease 
monotonically, in order to provide a smooth transition to the inflationary regime, where 
the amplitude is monotonically decreasing. Fig. 3. 

2. The bi-spectrum does not show a local maximum, see Fig. 6. As discussed in Sec. 4.2, 
during the inflationary regime the bi-spectrum's envelope decreases monotonically, 
while it increases towards the end of the curvature-dominated regime, indicating a 
maximum in the intermediate one. 

Both of these anomalies are artefacts caused by using the approximation (4.64) for (7{t), 
which underestimates the damping of oscillations for large k. In App. B, we revisit our 
computational strategy and explain these anomalies. 

We conclude that oscillations in the power-spectrum are indeed monotonically decreas- 
ing through the three regime, while the oscillation amplitude of the bi-spectrum has a local 
maximum in the intermediate regime. Further, our analytic results are trustworthy at the 
lower end of the k-domain in (4.76). 

4.3.1 Another Case Study and Observational Prospects 

Using our toy model V{(j)) = m^i?ii^/2 with cf) = 16 and = \/87r x 10~^, let us estimate 
the maximally attainable amplitude of the power-spectrum correction induced by an initial 
excitation a{ti) = ao. Since the amplitude is monotonically decreasing, we compute the 
correction at the lower boundary of the k-domain, 

with A = 10. Further, the strongest constraint on ctq is (3.9), 

^o«"^^ (4.78) 



with S = 0.1. For example, if rria = lOOy^P^/S we get ctq <C 0.015; in the plots of Fig. 5 and 
6, (To is close to this boundary. 

Requiring that at least one a-oscillation is completed in the intermediate regime, 27r/mo- = 
T < tf — ti, leads to 

— > ^ (4.79) 

^JQ arcsinh (;g) — arcsinh (^j 

PS 14.2, (4.80) 
where we use the time interval in (4.65). 
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We can approximate the amplitude of the power-spectrum correction in (4.73) by eval- 
uating it at A: = /cmin, since it decreases monotonically with /c, 




(4.81) 

We focused on the sine term (it carries the larger amplitude), took arcsinh(l/j4) ~ 1/^, 
ignored corrections of order 0{1/A), used the upper bound on do in the second to last step 
and took |p| ~ 3/2 in the last one. Thus, a signal generated close to the lower limit in 
(4.80), leads to a correction less then a quarter of a percent and for the largest possible mass, 
nia- ~ 1, we have 

-^<C 0.57. (4.82) 

Hence, in the most optimistic case with the largest possible fio allowed within our framework 
and a large mass, we might observe an oscillation amplitude of a few percent. 



5 Conclusion 

Motivated by ubiquitous metastable vacua in the multi- dimensional landscape of string the- 
ory, we consider two-field open inflation as a model of the early universe, as proposed in 
[17]. Changes in the power-spectrum due to super-curvature fluctuations as well as a direct 
detection of curvature are known tell-tale signs of open inflation. In this article, we propose a 
third complementary probe: oscillations in correlation functions caused by oscillations in the 
tunnelling field a after tunnelling, but before slow-roll inflation. The oscillation frequency of 
corrections to the power-spectrum as well as the bi-spectrum is set by the inverse of the scale 
factor, enabling (in principle) a direct detection of the background evolution and thus the 
curvature dominated regime. To use resonances induced by oscillating heavy fields as direct 
probes of the background evolution was first put forth by Chen [19], for instance, to discrim- 
inate infiation from proposed alternatives, such as bouncing cosmologies; since his results are 
not directly applicable to a curvature dominate universe, we revisit this computation in this 
article. 

We follow Chen's approach to compute correlation functions in three regimes: curva- 
ture domination, infiation and an intermediate regime smoothly connecting the former two. 
After specifying the smooth background, we compute oscillatory corrections to the Hubble 
parameter and its derivatives induced by oscillations in cr. At the perturbed level, we ignore 
fluctuations in a, paying special attention to the conditions under which such a truncation is 
feasible, leading to a strict upper bound on the allowed oscillation amplitude do in (3.9) and 
(3.10). This bound is not always imposed in Chen's work [19, 21], but satisfied throughout 
this article, leading to corresponding limits on signals in correlation functions. For larger ao, 
which are physically possible, fluctuations in all fields need to be included, since they are 
strongly coupled to each other. 

With this caveat in mind, we compute correlation functions via the in-in formalism; 
the resulting integrals allow for resonances between mode-functions and oscillations in back- 
ground quantities. Potentially resonating modes are deep inside the horizon (k ^> aH) where 
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differences between the Bunch Davies state and the (excited) state after the tunnehing event 
are suppressed; hence, we can take the former one as an approximation. We evaluate the 
resulting integrals via the stationary phase approximation, leading to analytic expressions 
for the power-spectrum correction and the bi-spectrum. The relevant expressions are given 
in (4.12), (4.47), (4.73) and in (4.14), (4.48), (4.75) for signals generated during the curva- 
ture, inflationary and intermediate regime respectively. We recover Chen's results [19] in the 
inflationary regime. 

We provide concrete case studies during curvature domination, Sec. 4.1.1, and the in- 
termediate regime, Sec. 4.3.1, to discuss the amplitude of signals and the feasibility of an 
observation. We find that signals caused by the initial excitation of a directly after the tun- 
nelling event are too weak (within our computational strategy) to be detectable with current 
experiments. Larger signals caused by a bigger initial amplitude of a after the tunnelling 
event are conceivable, but require a full multi-field treatment at the perturbed level. How- 
ever, a subsequent excitation can generate oscillations on top of the power-spectrum of a 
few percent, see (4.28) and (4.82), without violating the applicability of our approximation 
scheme. To achieve such large signals, the initial amplitude of a as well as its mass need 
to be relatively large, and the signal needs to be generated shortly before or during the 
intermediate regime. Furthermore, slow-roll inflation should not last much longer than the 
minimal amount needed to bring down spatial curvature below observational limits. If several 
oscillations on top of the power-spectrum were detected, the distinct change in the frequency 
uj{k) could act as a fingerprint of open inflation. Regarding the detectability of oscillations 
on top of the CMB power-spectrum by WMAP and PLANCK we refer the interested reader 
to recent articles such as [21, 75-77]. 

In discussing the intermediate regime, we found a spurious local maximum in the power- 
spectrum, which turned out to be an artefact caused by our analytic approximations. A 
complementary direct numerical integration leads to a continuously decreasing amplitude of 
the oscillatory signal. On the other hand, a local maximum is present in the shape function, 
which pinpoints the transition region from curvature domination to inflation. However, the 
amplitude of the bi-spectrum is suppressed. The latter may be increased in models that 
include a direct coupling between fields. 

To conclude, resonances caused by oscillating, heavy fields lead to oscillations on top 
of the power-spectrum and a distinct bi-spectrum, which can indeed provide a complemen- 
tary, direct probe of two-field open infiation, if they were observed. A detection of open 
infiation would offer corroborating evidence for the string landscape and thus a window of 
opportunity to test one of the few generic predictions of string theory, see e.g. [7]. Since the 
scientific gain of a detection is large, it appears prudent to us to intensify searches for such 
signals in the data, in conjunction with searches for super-curvature fluctuations as well as 
spatial curvature, even if the conditions under which such signals are detectable with current 
experiments are not the most generic ones. 
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A Normalisation of the shape function 

We briefly discuss the normalization of the shape function, since different choices can be 
found in the literature. Recall the definition of S in (3.16), 

(C=^> = S{h, k2, A;3)^(27r)7<53(ki + ks + kg) , (A.l) 

i 

and note that the normalisation factor of has its origin in the definition of the two-point 
correlation function {^(^^) in (3.14), 

(C'> = ^(2vr)^<5(ki + k2). (A.2) 

Hence, the proportionality to (A;2fc2^3)^ of the shape function in (3.35), 

c 



S (fci, /C2, /ca) = (27r)4p2 ^ j ^^'^ ^(e^)^^fci'"fc2 + ^ Pe™. + c.c. , (A.3) 



is determined by the choice of normalisation and has no intrinsic meaning. For example, if 
we chose to normalise S with (k) / {kik2k3) instead of P'^{k) / {kik2k^)'^ , we end up with a 
shape function, usually called of the form 



MkiMM) = j:^^i / dTG3A(e?7X4^-^ + 2 perm. + C.C., (A.4) 

V / C ''''"min 



c 

which simplifies for a Bunch-Davies mode function to 



A{K) = I ^ ^ {-2Kcos{Kr)+6f{T)sm{KT)) dr. (A.5) 

B AnomaUes in the Analytic Approximations 

We wish to elaborate on the origin of the two unexpected features in our analytic results: the 
presence of a local minimum in the oscillation-amplitude of the power-spectrum correction 
and the absence of a local maximum in the oscillation-amplitude of the bi-spectrum in the 
intermediate regime. 

B.l Power-Spectrum Correction 

Since we performed an expansion around the time of equality and the threshold wave number 
is A; ^ keq, it is conceivable that the observed increase is merely an artefact. To test this 
hypothesis, we complement our analytic estimates by direct numerical integration. 

Firstly, we compute numerically the solution to the full equation of motion (2.5) for a 
with ifbgr = a/p0/3 coth(t y^p^/S) . To separate -ffo-osci into smooth and oscillatory compo- 
nents, we take a moving average over a sufficient number of time steps, covering more than 
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Figure 7: Numerical computation of the power-spectrum correction with = 75 \J p^/'i 
and (To = 0.01. A monotonicahy decreasing amplitude of the oscillations is visible (compare 
to Fig. 5), leading to a smooth transition to the inflationary regime, Fig. 3. However, the 
numerical result is contaminated by our imperfect subtraction of the smooth contribution, 
leading to irregularities in the oscillations. 



one oscillation. We obtain the slow-roll parameters, substitute them into (3.34) and integrate 
numerically to compute the correction to the power-spectrum, see Fig. 7 for a representative 
example. 

While the result is contaminated by our imperfect split into smooth and oscillatory 
components, we clearly observe a monotonically decreasing amplitude, confirming our suspi- 
cion that the minimum in the amplitude is indeed an artefact of the approximations used. 
To see exactly where the problem arises, consider the integral (3.34). The key element is the 
contribution Ae/ebgr, which simplifies to 



^bgr ^^bgr 

where i/osci is the only unknown; according to the approximation in (4.64) as well as the 
numerical results, this function can be decomposed into a small oscillatory component with 
frequency a;(t), overlaid onto a smoothly decaying envelope A{t), 

Ae/ebgr = A(t)/iirbgr sm{uj{t)) , (B.2) 

as shown in Fig. 8 (without loss of generality we can choose A{t) > 0). Thus, the power- 
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Figure 8: The oscillating component of the slow-roll parameter, computed numerically for 
TTio- = 100 and ^ p^j'i. The deviation from Ae/ebgr = of the average is due to our imperfect 
subtraction of the leading order smooth component, which has no impact on our conclusions. 



spectrum correction becomes 



10 



— (fc)~ / (sin(2A;r -a;(t)) -hsin(2A;r-ht^(t)))di. (B.3) 

f ■ 

Here, the extra factor of "a" in the denominator stems from the transformation to physical 
time. Due to the restrictions on the r-domain (r < 0), only one of the sine terms contributes 
to the asymptotic behaviour of the integral. Since a minus sign can always be reabsorbed 
into the definition of uj{t) by exploiting the periodicity of the sine- function, we can freely 
take the contributing term to be the first one. Using the method of stationary phase, the 
envelope of the resulting oscillations in the power-spectrum is given by 



with 



(B.5) 



and the solution of /c* = a{duj/dt)/2 = g ^{t^f) = g^"^ . Substitution of the background 
variables yields 

AP 

-^(^*)cnv ~ KA{g^){-a{g*))x{g*y^^'^ ■ (B.6) 

Let's examine this expression carefully. Intuitively, the envelope should be a mono- 
tonically decreasing function of A;*, so as to smoothly connect with the beginning of the 
inflationary regime. Since monotonicity is preserved by composition but not by multiplica- 
tion of functions, we note the origin and behaviour of the different multiplicative factors: 

• k-t comes from the mode function, i.e. the Bunch-Davies vacuum oscillations; it trivially 
increases with /c*. 



For simplicity, we neglect the contributions of the /-term in the power-spectrum integral, (3.34). 
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Figure 9: rria = 50 ^ p^j?). The analytic approximation used in Sec. 4.3 begins to devi- 
ate from the fuh numerical result at late times (high wave numbers) exemplified by under- 
damping and a phase shift in the cr-oscillations, leading to the unexpected increase in the 
power-spectrum signal, see Fig. 5. 

• ~0'ig*) comes from the background interaction of the inflaton with the spatial curva- 
ture. It is monotonically increasing as a function of t*, so that its overall behaviour is 
determined by . 

• ^((7*) - the amplitude of the a field oscillations is naturally decreasing with physical 
time due to the damping effect of the Hubble term. Its behaviour in momentum space 
is again determined by g*. 

• x{d*)~^^'^ ' this term captures the strength of the resonance between the vacuum and 
cj-oscillations and its behaviour cannot be determined without a specific model. 

To summarise, the monotonicity of the power-spectrum correction is determined primar- 
ily by the behaviour of the resonating wave mode k^: = g^^ and the resonance mechanism x*- 
In the case of the approximation (4.64), both g^, and x(S'*)~^^^ increase monotonically. So, the 
only decreasing factor in (B.6) is the one related to the oscillation amplitude, A{g^). We con- 
clude that the break of monotonicity observed in Fig. 5 has to be due to the insufficient damp- 
ing of the cr-oscillations at high wave numbers (late times) in the power-law model. Indeed, 
(4.63) can be asymptotically expanded to reveal a factor of Exp (f (\/2 -|- {tf^ — t) y^p^T^)^) , 
which enhances the damping at late times. This is illustrated in Fig. 9, which shows a 
comparison of the fi-oscillations obtained numerically and our analytic approximation. 

We conclude that the power-spectrum correction does indeed monotonically decrease in 
the intermediate regime; only the approximation (4.73) breaks down beyond the stationary 
point of its envelope, leading to the spurious increase visible in Fig. 5. 

B.2 Bi-spectrum 

Applying the argument outlined above to the bi-spectrum dampens the monotonic increase of 
the amplitude at higher wave numbers as well, conceivably leading to the expected maximum. 
As the position and height of this maximum depends sensitively on the chosen shape of the 
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slow-roll potential, we are satisfied with our current understanding and leave a numerical 
computation to the future. 
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